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Abstract 

Gortler vortices arise in boundary layers along concave surfaces due to 
centrifugal effects. This paper presents some results of an experiment 
conducted to study the development of these vortices on an airfoil with 
a pressure gradient in the concave region where an attached laminar 
boundary layer was insured with suction through a perforated panel. 
The sublimating chemical technique was used to visualize Gortler 
vortices and the velocity field was measured by laser velocimetry. The 
vortex wavelength clearly varied with Gortler number as predicted by 
linear theory. Both flow visualization and velocity measurements 
indicated vortex damping in the convex zone. Secondary instability 
was observed at the higher Gortler numbers. 
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Symbols 


c 

p 

pressure coefficient 

c 

chord length 

G US 'F 

v I r 

Gortler number 

g m 

Gortler number at maximum curvature 

^00 

free-stream Mach number 


Reynolds number based on chord 

« u - 
v 

unit Reynolds number 

r 

radius of curvature 


free-stream velocity 

u,v,w 

velocity components in laser coordinates 

u 

streamwise boundary-layer velocity component 

x, y, z 

model coordinates 

a 

wave number 

5 P 

r \U 

boundary layer thickness parameter 

A 

V 

nondimensional wavelength 


Introduction 


Gortler vortices arise in boundary layers along concave surfaces due to 
centrifugal effects. These counter-rotating streamwise vortices are 
one of three known flow instabilities which lead to boundary layer 
transition. Coupled with Tollmien-Schlichting waves and crossflow 
vortices these Gortler vortices can trigger early transition to 
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turbulence. Gortler vortices may play an important role in internal 
flows with concave curvature and on modern supercritical laminar- 
flow-control wings which have concave regions near the leading and 
trailing edges of the lower surfaces. 

This centrifugal instability was first treated analytically by Gortier 1 
after whom the vortices are named. The stability is governed by the 
parameter 



where 8 r is a suitable boundary layer thickness, r is the radius of 
curvature, and v is the fluid kinematic viscosity. Many studies have 
since been devoted to the improvement and extension of Gortler’s 
analysis and a thorough survey of these efforts has been reported by 
Herbert. 2 There are a number of fundamental differences in these 
theoretical studies as to the details of the formulation of the problem as 
well as in the computed results. These are discussed in Floryan’s 
thesis. 3 The influence of external disturbances is analyzed by Floryan 
and Saric 4 and the influence of pressure gradient is discussed by Hall. 5 

There have also been a number of experimental studies of this problem, 
mostly conducted in curved channel flows. Tani et aZ, 6 " 9 attempted to 
correlate the streamwise vortex wavelength with wall curvature and 
free stream velocity (i.e., Gortler number) as predicted by theory. The 
failure to find such a correlation led them to the conclusion that the 
vortex wavelength is determined more by the particular edge effects of 
the experimental apparatus and less by the curvature of the model or 
the velocity of the oncoming flow. Several other experimenters 
including Bippes, McCormack et al, Winto et al, Yurchenco 
et aZ, 13 Babenko 14 and Yurchenco 15 , and Swearingen and 
Blackwelder 16 have also found that the wavelength is essentially 
independent of the Gortler number. Bippes 10 did find that when the 
stream turbulence was made isotropic, the most amplified disturbance 
according to theory could be generated. Experimental data from these 
studies seem to indicate that some preferred wave number disturbance 
is selected initially from those unstable ones and is amplified with 
increasing Gortler number while roughly following lines of constant 

U r X 1 

A 

v r 

where X is the wavelength. The present paper evaluates this parameter 
on the basis of tests conducted on a two-dimensional airfoil model. 
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Results from flow visualization tests using sublimating chemicals and 
laser velocimeter measurements are compared. 


Test Apparatus 

A 1.83-meter chord airfoil model was tested in the NASA Langley Low- 
Turbulence Pressure Tunnel (LTPT). A photograph of the model is 
shown in figure 1 . The model consists of two parts a structural 
element and a test element. The test element includes the leading edge 
and the upper surface back to midchord. The structural element 
consists of the spar and the remainder of the airfoil surface including a 
10% chord flap. With this design, several test region geometries can be 
examined. The present test element includes a concave region from 
x/c = 0.175 to x/c = 0.275. The minimum radius of curvature is 0.24 m. 
Attached laminar boundary-layer flow is insured by means of suction 
through a 0.11 x 0.76-m perforated titanium panel located in the 
compression part of the concave region. The suction region is divided 
into three spanwise suction strips. The suction in each strip is 
independently controlled by its own needle valve. The model surface is 
painted flat black over gel coat and fiberglass. Figure 2 shows the 
model pressure distribution. The 10% chord flap is used to control the 
stagnation point location to maintain a flat or slightly favorable 
pressure gradient ahead of the concave region. 

The LTPT is a pressurized, closed-circuit, continuous flow wind 
tunnel. The test section is 2.29 m high, 2.29 m long, and 0.91m wide. 
The contraction ratio is 17.6. The tunnel has excellent flow quality 
due, in part, to the nine screens in the settling chamber. The velocity 
fluctuations in the test section were found to be 0.025 percent at 0.05 
Mach number. 1 The pressure fluctuations at the test section wall, 
normalized with respect to free-stream pressure were 10" 5 at this 
Mach number. 

The present experiment was conducted at atmospheric pressure. The 
chord Reynolds number was varied from 1.0 million to 5.9 million and 
the Mach number ranged from 0.024to 0.125 yielding a Gortler number 
of up to 46. 

A specialized single axis three-component laser velocimeter was used 
to study the flow field in the test region (fig. 3). The single axis, five- 
beam optical configuration 1 uses the standard two color, two- 
component beam pattern with the two green beams (514.5 nm) 
arranged in the horizontal plane and the two blue beams (488.0 nm) 
arranged in the vertical plane forming a diamond pattern to measure 
the U or streamwise component (green beams) and the V or vertical 
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component (blue beams). At hird green beam is placed along the optical 
axis bisecting the angle between the original two green beams. The 
addition of this beam creates two additional fringe patterns within the 
sample volume. These additional fringe patterns are inclined 
symmetrically about the optical axis yielding equal contributions of 
the U component and equal but opposite contributions of the W or 
spanwise velocity component. The W velocity component is then 
obtained from the difference between the two signal frequencies. This 
is achieved by incorporating Bragg cells in the two outside beams to 
separate by frequency the three signals obtained from the three green 
fringe patterns. A filter network isolates the U component signal 
frequency and the two signal frequencies from the inclined fringe 
patterns which are then input to an electronic double balanced mixer. 
The lower frequency signal from the mixer is the frequency difference 
between the two input signals and is the W component signal frequency 
with a known bias depending on Bragg cell frequencies. The flow was 
seeded with tridecane particles using a particle generator located 
upstream of the settling chamber screens. The LV data were recorded 
and processed by a dedicated computer system. The velocity 
components were measured in the laser velocimeter coordinate system 
which looked down on the model at an angle of 11 . These velocity 
components were then transformed into wind-tunnel and model 
coordinate systems for analysis. 


Results and Discussion 

A thin layer of solid white biphenyl material was sprayed over the black 
model surface to visualize the flow. This hydrocarbon material 
sublimes due to the heat transfer at the airfoil surface which is 
proportional to the surface shear stress 19 . The flow pattern is made 
visible due to the differential surface shear stress distribution under 
the layer of counter-rotating Gortler vortices. Elapsed times of about 
30 minutes to 1 hour were required for the pattern to emerge clearly, 
depending on the free-stream velocity. Representative flow patterns 
are shown in Figures 4 (a), (b), and (c) taken at chord Reynolds numbers 
of 2.24, 3.21, and 3.67 million respectively. The dark bands represent 
the high shear stress regions where the chemical layer has sublimated 
revealing the black background of the airfoil surface, whereas, the 
white bands correspond to low shear regions. A set of black and white 
bands constitutes a pair of counter-rotating vortices and represents 
the wavelength of these vortices. 

The photographs in figure 4 show the region 0.20 < x/c < 0.40 near the 
center of the model span. The perforated titanium suction panel is 
visible at the upstream edge of each photograph. No streamwise 
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streaks were observed in the region ahead of the concave zone. The 
streaks observed downstream of the suction panel are essentially 
uniformly spaced along the span. Wavelength was determined by 
taking the average of the number of pairs of streaks over a 15 to 45 cm 
span. The nondimensional wavelength parameter A was computed 
using this value of dimensional wavelength and the mean flow and 
geometric parameters. The dimensional wavelength appears to remain 
constant throughout the concave region. Beyond the concave region 
the streaks decrease considerably in contrast but remain visible back to 
the jagged transition line. This indicates damping of the Gortler 
vortices in the convex region which is confirmed by the laser 
velocimeter measurements discussed below. Furthermore, in the 
convex zone, occasional vortex mergers are visible. 

Figures 5 and 6 show the variation of the streamwise velocity 
component along the span from 1.5 to 4.5 cm left of midspan at several 
heights above the model surface for chord Reynolds numbers of 1.0 and 
2.1 million, respectively. These Reynolds numbers correspond to Mach 
numbers of 0.024 and 0.05. The measurements presented are for chord 
locations of x/c = 0.15, 0.175, 0.20, 0.225, 0.25, 0.275, 0.30, and 0.3375 
for M x - 0.05 and all of these except x/c = 0.15 for M x — 0.024. The 
first and last of these locations are just ahead and aft of the concave 
zone, respectively. All three boundary-layer velocity components were 
measured, but only the streamwise boundary-layer component data are 
shown since this component is much larger than the other two. The 
spanwise scan lines for the higher Reynolds number data were taken in 
a vertical plane. The scanning method was later changed so that the 
scan lines for the lower Reynolds number data lie in a plane normal to 
the model surface. Figures 5 and 6 show a general decrease in velocity 
with increasing 2 . This is caused by minute deflections of the floor 
supporting the heavy laser optics traversing mechanism. The heights y 
indicated on the figure keys are correct at 2 = 1.5 cm. The maximum 
deflection of the laser control volume was measured to be 
Ay = 0.018 cm. The concave region of the model begins at x/c = 0.175. 
However, figures 5(a), 6(a), and 6(b) show small, but definite, periodic 
velocity variations along the span at x/c = 0.175. The disturbance 
wavelength was determined by direct examination of the streamwise 
velocity plots and by examination of the autocovariance function 
determined using the Fast Fourier Transform (FFT) technique. 

For G m - 29.9, figure 5(d) shows that the periodic velocity pattern is 
well developed by x/c = 0.25. The increase of the disturbance amplitude 
from x/c = 0.175 to 0.275 is clearly evident (figs. 5(a) - (e)). In the 
convex zone beyond x/c = 0.275, the disturbance amplitude decreases 
(fig. 5(f) and 5(g)). The dimensional disturbance wavelength remains 
essentially constant at X = 0.6 cm throughout the measurement range 
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at all heights except possibly near the boundary-layer edge at 
x/c = 0.275 where a doubling or tripling of the wavelength may be 
beginning. 

For G m = 36.0, the disturbance velocity pattern is more complex. 
Periodic spanwise variations are clearly evident from x/c = 0.15, 
(fig. 6(a)) and are well developed at the lower levels of the boundary 
layer by x/c = 0.225 (fig. 6(d)). However, a strong secondary 
disturbance with wavelength 4 to 5 times that of the basic instability 
abruptly appears at x/c = 0.225 (fig. 6(d)). The maximum velocity 
difference reaches sixty percent of the local edge velocity for y = 0.2 cm 
(fig. 6(d)). This secondary instability has been noted before by 
Wortman 20 and Bippes 10 . The secondary disturbance decays markedly 
in the accelerated flow region beyond x/c = 0.225 and the following 
convex region (fig. 6(d)-6(h)). Visual examination of the velocity plots 
and the autocovariance functions for the data of figure 6 indicate a 
disturbance wavelength of X - 0.3 cm. This wavelength is clearly 
evident from x/c = 0.225 to 0.30. At x/c = 0.3375, the dominant 
wavelength appears to be X = 0.6 cm or twice the basic wavelength; 
however, the autocovariance functions for this data indicate that 
variations of the shorter wavelength are still present. 

Laser velocimeter measurements of streamwise velocity at all chord 
locations as well as at all heights above the surface were used to 
determine disturbance wavelength, using autocovariance function. 
Except where secondary instability was present, the wavelength did 
not vary bymorethan 10% throughout the measurement zone at a given 
free stream condition. Wavelengths, determined from laser 
measurements, agree well with those observed in flow visualization. As 
in all previous experiments, the dimensional wavelength of the 
disturbance was conserved in the flow direction but unlike earlier 
experiments the wavelength was observed to vary appreciably with free 
stream conditions. The variation in wavelength was most noticeable at 
lower speeds as seen from figure 7 where X is plotted against Gortler 
number. The theoretical curve corresponds to wavelength computed 
from maximum amplification conditions. 3 It is seen that experimental 
values of X lie very close to wavelengths predicted on the basis of linear 
theory. The nondimensional wavelength parameter A obtained from 
flow visualization and laser velocimeter measurements is also shown 
plotted as a function of Gortler number G M in figure 8 where it is 
compared with existing experimental data and Floryan’s 3 theoretical 
results. 
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Conclusions 


1. A fixed, essentially uniform, vortex spacing was observed in the 
concave zone by both flow visualization and laser velocimeter 
measurements for each flow condition. 

2. Both flow visualization and laser velocimeter measurements 
show a vortex wavelength which varied with Gortler number in 
accordance with linear stability theory. 

3. A significant, abrupt decrease in streak contrast indicated vortex 
damping in the convex zone. 

4. The velocity measurements showed both disturbance 
amplification in the concave zone and damping in the following 
convex zone. 

5. Flow visualization did not indicate a vortex pattern ahead of the 
concave zone; however, velocity measurements indicated small 
spanwise variations of the streamwise velocity component at or 
just ahead of the beginning of the concave zone. The wavelength 
of this variation matched the observed wavelength in the 
downstream concave region. 

6. The flow visualization showed some apparent vortex mergers in 
the convex zone. 

7. The velocity measurements indicated a stationary secondary 
instability in the concave zone. 
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Figure 1.- Gortler model in LTPT. 




Figure 2.- Model schematic diagram and pressure distribution. 
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Figure 3.- Laser velocimeter optics 




2 


I 


— 




Figure 4(a).- Flow visualization using sublimating chemicals M = 0.05, 
R c = 2.24xl0 6 . 
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Figure 4(b).- Flow visualization using sublimating chemicals M = 0.075, 
R c = 3 .2 1 x 1 0 6 . 



Figure 4(c).- Flow visualization using sublimating chemicals M = 0.089, 
R c = 3.67xl0 6 . 
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Figure 5(a).- Spanwise variation of streamwise velocity component 
R c = 1 .0x10 , M x = 0.024, G = 29.9, x/c = 0.175. 



Figure 5(b).- Spanwise variation of streamwise velocity component 
R c = 1.0x10 , M x = 0.024, G = 29.9, x/c = 0.20. 
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Figure 5(c).- Spanwise variation of streamwise velocity component 
R c = 1.0x10®, M x = 0.024, G = 29.9, x/c = 0.225. 



Figure 5(d).- Spanwise variation of streamwise velocity component 
R c = 1.0x10 , M x = 0.024, G = 29.9, x/c = 0.25. 
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Figure 5(e).- Spanwise variation of streamwise velocity component 
R c = 1 .0x10 , M x = 0.024, G = 29.9, x/c = 0.275. 



Figure 5(f).- Spanwise variation of streamwise velocity component 
R c = 1 .0x10 , M x = 0.024, G = 29.9, x/c = 0.30. 


16 






Figure 5(g).- Spanwise variation of streamwise velocity component 
R c = 1 .0x10 , M x = 0.024, G = 29.9, x/c = 0.3375. 



Figure 6(a).- Spanwise variation of streamwise velocity component 
R c = 2.1x10®, M oo = 0.05, G =36.0, x/c = 0.15. 
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Figure 6(b).- Spanwise variation of streamwise velocity component 
R c = 2.1x10 , M o0 = 0.05, G =36.0, jc/c = 0.175. 



Figure 6(c).- Spanwise variation of streamwise velocity component 
R c = 2.1x10®, M o0 = 0.05, G =36.0, x/c = 0.20. 
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Figure 6(d).- Spanwise variation of streamwise velocity component 
R c = 2.1x10 , M oo = 0.05, G =36.0, x/c = 0.225. 



Figure 6(e).- Spanwise variation of streamwise velocity component 
R c = 2.1x10®, M oo = 0.05, G =36.0, x/c = 0.25. 
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Figure 6(f).- Spanwise variation of streamwise velocity component 
R c = 2. 1x10°, M oo = 0.05, G =36.0, x/c = 0.275. 



Figure 6(g).- Spanwise variation of streamwise velocity component 
R c = 2Ax\0 6 , M oo = 0.05, G =36.0, x/c = 0.30. 
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Figure 6(h).- Spanwise variation of streamwise velocity component 
R c = 2.1x10 , M oo = 0.05, G =36.0, x/c = 0.3375. 



Figure 7.- Variation of dimensional wavelength with Gortler number, 
comparison with theory. 
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Figure 8.- Variation of nondimensional wavelength parameter with Gortler 
number, comparison with theory and other experiments. 
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